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.ABSTRACT 

,  C 

'“"'This  thesis  presents  an  analysis  of  stochastic  duels  in- 

/ 

volving  two  opposing  weapon  systems  with  constant  rates  of 
fire.  The  duel  was  developed  as  a  stationary  Markov  chain 
with  stochastic  matrices  of  transition  probabilities  con¬ 
structed  from  the  single  shot  kill  probabilities  of  the 
weapon  systems.  A  comparison  was  made  of  the  presented 
Markov  chain  analysis  results  with  results,  from  other  ac¬ 
cepted  conditional  probability  methods.  As  expected,  this 
comparison  established  the  validity  of  the  Markov  chain  anal¬ 
ysis  and  indicated  advantages  of  the  Markov  chain  approach  in 
analysis  of  discrete  process  stochastic  duels.  The  analysis 
was  then  extended  to  the  two  versus  one  duel  where  the  three 
weapon  systems  were  assumed  to  have  fixed  rates  of  fire.  ^  j 

t 
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I INTRODUCTION 


A.  BACKGROUND 

The  theory  of  stochartic  duels  has  been  developed  in 
order  to  evaluate  the  probability  of  a  weapon  system  sur¬ 
viving  an  engagement  with  anothe..  weapon  system.  The  effect 
of  weapon  system  parameters ,  such  as  single  shot  kill  proba¬ 
bilities,  rate  of  fire,  and  tactual  (time)  advantage  can 
be  determined  in  the  analysis  of  these  duels. 

Joseph  J.  Schoderbeck  (1962)  and  Trevor  Williams  and 
C.  J.  Anckar,  Jr.  (1963)  developed  the  basic  theory  of  the 
stochastic  duel  and  the  "fundamental  duel"  was  defined  as 
follows : 

1.  Two  comb  .tants,  A  and  B,  fixed  at  each  other  until 
one  was  killed. 

2.  The  time  between  rounds  fired  was  either  a  constant 
or  a  random  variable  of  known  but  different  density  function 
for  each  combatant. 

3.  Each  combatant  had  a  different'  known  but  fixed  single 
shot  kill  probability. 

4.  The  duel  began  with  each  combatant  having  unlimited 
ammunition  supplies. 

5.  Both  combatants  had  unlimited  time  in  which  to  score 
a  kill. 

6.  A  tactical  (time)  advantage  in  firing  the  first  shot 
was  assigned  to  one  of  the  combatants. 
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Since  the  first  development  of  stochastic  duel  theory  in 
1962,  many  extensions  to  the  theory  have  evolved.  Analysis 
of  weapon  system  duels  constrained  by  ammunition  limits,  time 
limits,  or  varying  single  shot  kill  probabilities  has  been 
accomplished.  Also,  distributions  of  the  time  to  kill  and  of 
the  number  of  rounds  fxred  have  been  determined.  In  addition 
two-versus-one  duels  and  two-versus-two  duels  have  been  de¬ 
veloped.  As  more  combatants  participate  in  the  duel,  it  was 
shown  that  the  results  were  in  keeping  with  Lanchester  models. 

The  objective  of  this  thesis  was  to  analyze  the  "funda¬ 
mental  duel"  involving  fixed  time  between  firings,  using 
stationary  Markov  chains.  As  expected  the  Markov  process 
analysis  was  shown  to  be  equivalent  to  the  conditional  proba¬ 
bility  methods.  But  in  the  Markov  analysis  of  the  stochastic 
duel,  there  existed  advantages  over  other  analysis  methods 
that  are  indicated  in  Section  II.  Finally,  the  Markov  anal¬ 
ysis  was  extended  to  the  two-versus-one  duel  where  all  three 
weapon  systems  were  assumed  to  have  fixed  rates  of  fire. 

B,  PREVIOUSLY  DEVELOPED  ANALYSIS 

Stochastic  duels  have  been  analyzed  utilizing  conditional 
probability  methods.  Schoderbeck  [Ref.  1]  was  first  respons¬ 
ible  for  the  analysis  of  the  fundamental  duel  involving  fixed 
time  between  firings.  A  summary  of  his  analysis  is  presented 
below. 

A  (friendly  force)  and  B  (enemy  force)  each  possessed  a 
single  weapon  system  with  single  shot  kill  probabilities  PA 
and  pB  respectively. 
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Let  a  =  time  between  A's  firings, 
b  =  time  between  B's  firings. 

It  was.  assumed  that  at  time,  t-0,  A  fired  his  first  shot 
at  B  and  that  sometime  1  r,  T,  B  returned  fire.  Also,  it 
was  assumed  that  a=b,  that  T^a,  and  that  each  firing  was  an 
independent  event.  Then  the  firing  sequence  looked  like  that 
depicted  in  Fig.  1. 


A 


B 


«-a-' 

— X - X - X - X - X - 

time 

«-b+ 

— - - X - X - X - X - 

t=0 

Figure  1.  Firing  Sequence  of  A  and  B. 


Let  Pr[A(tQ)]  =  Prob[A  was  alive  and  B  was  dead 
at  time  t=tQ]  and 

Pr[A(«)]  =  ^~PrtA(ton* 

Considered  first  was  t,\-*  interval  a<tQ<a+T.  Within  this 
interval  there  was  the  possibility  that 

a)  A  fired  exactly  two  times. 

b)  B  fired  exactly  one  time. 

Thus  A  was  alive  and  B  was  dead  if  and  only  if 

a)  B  was  killed  on  A's  first  'shot  (so  did  not  fire 
a  shot) ,  or 

b)  B  missed  with  his  one  shot  and  was  killed  by  A's 
second  shot. 

But,  Pr[B  was  killed  on  A's  first  shot]  =  pa 

Pr[B  missed  with  his  first  shot  and.  was  killed  with 


A's  second  shot] 


<1-PA)  U-PB)PA 


Wa 


and  thus,  t>r[A(to)]  =  PA+(IA<1BPA  for  a<tQ<a  +  T* 

Then  considered  was  the  interval  a+T<t  <2a+T.  In  this 

o 

interval  there  was  the  possibility  that 

a)  A  fired  exactly  three  shots . 

b)  B  fired  exactly  two  shots. 

And  here  A  was  alive  and  B  was  dead  if  and  only  if 

a)  B  did  not  fire  at  all  (killed  by  A's  first  shot) , 
or  b)  B  fired  once,  missed,  and  was  killed  by  A's  second 
shot,  or 

c)  B  fired  twice,  missed  twice,  and  was  killed  by  A's 
third  shot. 

But,  Pr[B  did  not  fire  at  all]  =  pA 

Pr[B  fired  once,  missed,  and  was  killed  by  A's  second  shot] 

-  (1-PA}  (1'VPa  =  qAqBPA 

Pr[B  fired  twice,  missed  twice,  and  was  killed  by  A's 
third  shot]  =  (1“PA)  (1“PB)  (1“PA)  d-PB)?A  =  <3AqBPA 

Hence,  Pr[A(tQ)]  =  pA+qAqBpA+qA2qB2pA  for  a+T<tQ<2a+T 
By  induction  then  for  the  interval  (n-1) a+T<tQ<na+T 
PrtA(t0>]  -  PA+PA<9a9B)+t,A(qAqB)2+---+PA(qA<3Bin 
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Since  p^  >  0  and  PB  >  0  implied  qA<l  and  qB<l  and 


lim 

n->oo 


<qAqB) 


n 


0 


Ancker  and  Williams  [Ref.  2]  developed  a  method  for  anal¬ 
yzing  the  fundamental  duel  with  fixed  firing  rates,  where 
a/b  but  T=0.  A  summary  of  their  analysis  is  presented  below. 

The  ratio  jj  was  assumed  to  be  a  rational  number  and  if 
a  and  b  contained  a  common  factor,  this  ratio  was  reduced  to 
a/3  where  a  and  3  were  relatively  prime  integers. 

Let  n  =  the  largest  number  of  times  3  was  contained 

in  a. 

Let  r  =  the  remainder  when  3  was  divided  into  a,  then 


a=n3+r. 


Let  P (A)  =  Pr(A  won  the  duel)  =  Pr(A  was  alive,  B  was 
dead) .  Then  A  won  on  the  jth  shot  if  he  missed  B  on  his  first 
j-1- shots  and  had  hit  on  his  j  round,  while  B  missed  with 
his  first  k  rounds  where  k=[j(a/3)l  and 

[X]  *  largest  integer  less  than  or  equal  to  X. 

Then  the  probability  that  A  won  the  duel  was 


P<A)  "  jh* AV  “qB 


PA  9b 


* 

j-l„  jn+[ j (r/3) 1 


PA  jc 

i  3*  a  • 
iXqB  3=0 


[(j+l)  a/p] 
B 


(1-3) 
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C.  APPLICATIONS  OF  STOCHASTIC  DUELS 

As  analytical  formulations  of  combat  operations, i stochas¬ 
tic  duels  are  of  great  value  in  providing  insight  for  the 
design  of  new  weapon  systems  in  advanced  time  frames . 
Stochastic  duels  lend  themselves  easily  to  parametric  anal¬ 
ysis.  By  varying  parameters  such  as  single  shot  kill  proba¬ 
bilities  or  .firing  rates,  it  is  possible  using  this  analysis 
to  determine  tradeoffs  between  volume  of  fire  versus  accuracy 
of  fire. 

The  military  can  apply  the  analysis  of  stochastic  duels 
to  missile  sites,  artillery,  torpedoes,  and  most  other  weapon 
systems.  Using  realistic  parameters,  the  survival  proba¬ 
bilities  can  be  determined  in  engagements  with  enemy  weapon 
systems . 

Analysis  of  duels  involving  fixed  times  between  firings 

can  also  be  compared  with  analysis  of  duels  involving  random. 

times  between  firings,  using  the  same  parameter  values.  This 

comparison,  could  yield  useful  information  concerning  the 

circumstances  where  a  fixed  rate  of  fire  was  more  effective 
✓ 

than  a  variable  rate  of  fire. 

Section  II  contains  the  presentation  of  the  stationary 
Markov  chain  analysis  of  the  fixed  rate  of  fire  stochastic 
duel.  This  method  was  verified  also  in  Section  II  by  com¬ 
paring  results  with  conditional  probability  analysis  results 
which  were  presented  in  this  Section.  Examples  of  Markov 
analysis  of  some  duels  completed  Section  II. 
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Section  III  presented  an  extension  of  the  analysis  to 
the  two-versus-one  duel,  including  an  example. 

The  conclusions  reached  after  the  analysis,  and  recom¬ 
mendations  for  further  study  are  contained  in  Section  IV. 
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11 •  STOCHASTIC  DUEL  -  FIXED  RATES  OF  FIRE 

A.  DEVELOPMENT  OF  THE  MARKOV  CHAIN  MODEL 

A  and  B  represented  two  combatants,  each  possessing  one 
weapon  system.  The  weapon  systems4  respective  single  shot 
kill,  probabilities  were  p^  and  pB.  Further,  it  was  assumed 

that  both  A  and  B  fired  at  a  fixed  rate  so  that  a  and  b 

represented  the  respective  times  between  firings  (in  seconds) 
of  A  and  B.  Then, 

\  60 

aA  =  —  =  A's  rate  of  fire  (shots/min.) 

a 

=  j-p-  =  B’s  rate  of  fire  (shots/min.) 

It  was  assumed  also  that  both  a  and  b  were  rational  numbers 
and  that  A  had  an  arbitrary  tactical  time  advantage,  T  (sec¬ 
onds)  .  Then: 

A  fired  his  first  shot  at  time,  t=0,  and 

B  fired  his  first  shot  at  time,  t=T. 

J.U 

Af  then,  fired  his  j*  shot  at  time 

tAj*  “  a  ( j*-l)  (II-l) 

and  B  fired  his  k^  shot  at  time 

tBk  -  T+b (k-1)  (II-2) 

illustrated  in  Fig.  2. 
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a  -> 


t=0 

Figure  2.  Time  Relation  Between  A  and  B's  Firings 


At  any  time  in  the  process,  the  duel  was  in  one  of  four 
possible  states: 

STATE  1:  A  was  alive;  R  was  alive 
STATE  2.i  ■  A  was  alive;  B  was  dead 
STATE  3:  A  was  dead;  B  was  alive 
STATE  4:  A  was  dead;  B  was  dead 

It  was  noted  that  STATE  4  occurred  only  if  A  and  B  fired 

simultaneously.  Whenever  A  or  B  (or  both)  fired,  the  duel 

underwent  a  transition  from  one  state  to  another.  Transition 

matrices  associated  with  the  three  possible  occurrences  that 

could  cause  a  transition  from  one  state  to  another  are  in  the 

form  of  the  following: 

1)  A  only  fired 


2)  B  only  fired 
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3)  A  and  B  fired  simultaneously 


qAqB  PAqB  qAPB  PAPB 


where  qA  =  l-pA  and  qB  =  l-pB  . 

The  initial  state  vector  at  the  beginning  of  the  duel  was 
P  *(1,0, 0,0)  since  the  duel  was  always  in  STATE  1  (A  and  B 
both  alive)  at  the  outset.  After  m  transitions,  however, 
the  state  vector  was  represented  by 

m  o  1  2  •  •  •  m 

where  P^,i=l,2, • • * ,m  was  PA,  Pfi,  or  PAB  depending  on  the  fir¬ 
ing  sequence.  This  firing  sequence  was  in  turn  dependent 
upon  T,  XA,  and  Xfi. 

To  compute  P^,  the  number  of  times  A  fired  in  the  period 
[0 ,T]  first  had  to  be  determined.  This  was  done  as  follows: 


Set  t. 


=  tT 


a(j*-l)  =  T  ,  and  solved  for  j* 

j*  s  —  +  1  and  if  [j*]  indicated  the  greatest 
Si 

T 

integer  less  them  j*  ,  then  j  ®  [j*]  *  (~  +  1]  repre¬ 
sented  the  number  of  times  A  fired  in  the  period  [0,T],  Then, 


if  •  P.  fr»r  rn<-] 

O  A  - 

P  ’P,3 *P . . i »P . . o •  »P«  for  j<m 

o  A  3+1  3+2  • • •  m  J 


(II-3) 


where  P 


it  i  =  j  +  1,  j  +  2,...,  m  was  PA,'Pfi  or  P; 
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depending  on  firing  sequence  derived  from  XA  and  Xfi.  Also 
beginning  with  the  sequences  of  matrices  were  periodic 

as  shown  below: 

A  fired  his  (j+l)st  shot  at  time,  t>T 
B  fired  his  1st  shot  at  time  t=T 
The  time  between  B's  1st  shot  and  A's  (j+l)st  shot  was 


bA,j+l  “  fcBl  ~  a3  “  T 


(II-4) 


If  t  represented  the  period  in  seconds  of  the  firing  sequence 
then  the  difference  in  time  between  the  shot  fired  by  A  at 
time  tA  +  nr  and  the  shot  fired  by  B  at  time 
tD,  +  m  ,  n=l,2,*.*,  equaled  the  time  interval  of  Eq.  II-4 . 
This  was  proven  by  the  following  method: 


—  and  Xn 
a2  B 


°1 

=  g—  where  a^,  a2,  b-^,  and  b2  were 
2 


integers  since  a  and  b  were  assumed  to  be  rational  numbers. 

Let  L  »  lowest  common  multiple  (L.C.M.)  of  a2  and  b2 

L  ~  L.C.M.  (d2 ,b2) 

Then  ~  and  were  both  integers. 

a2  b2 

LXA  =  ai(iJ} 

lxb  =  bl(5J) 

Let  L*  =  L.C.M.  (LXa,  LXfi) 


Then  t  =  v?~v*  (seconds)  was  the  period  of  the  firing 

lxaxb 

sequence.  A,  then,  fired  t(g^)  =  g  shots  in  t  seconds  and 

XB  T 

B  fired  t  (^)  =  g  shots  in  t  seconds. 
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Let  nA  =  number  of  the  shot  A  fired  at  time 

t.  ...  +  nx  and  nu  =  number  of  the  shot  B  fired  at  time 
A, 3+1  B 


tB1tnr  . 


nA  “  J  +  1  +  tn't<So'>1 


nQ  =  1  +  {nx  (—•)}  for  n  =  1,2, 


4-  Aa  60  . 

%nA  *  a(nA-l)  =  {j+nxC™)}  =  — 3-  +  nT 

A  A 


:B,n„  *  1  +  -  T  +  XT  <"*<«$>  =  T  +  nx 

3  3 


so  that  t_.  -  t.  „  -  r— 

B,nB  h'nh  AA 


60.  60. 

- — 1  +  nx  -  (T+nx)  =  r— -  T  =  a .  -  T 
aA  aA  j 


*  Vj+l  "  bBL 

In  a  cycle  then  the  total  number  of  shots  fired,  N,  was  rep¬ 
resented  by  N  =  x(g^)  +  !(£$•)  but  =  EX“"  and 

B 

*‘s§»  -  Q-  80  that  N  =  r-  (XJ  +  ^  • 

If  N*  represented  the  maximum  number  of  matrices  per  cycle, 


then, 


N*  = 


N,  (~J-  of  type  pA  ;  of  type  B)  if  there 
B  A 

was  no  simultaneous  firing 


^  N-l,  ((jjj.  -1)  of  type  PA,  -  1)  of  type  P0 

\  B  A 


and  1  of  type  PAfi)  if  there  was  a  simulta¬ 
neous  firing. 
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The  state  of  the  duel  after  n  transition  cycles  (instead 
of  m  transitions) ,  then,  was  written  as: 


P.  =  P . •  P  J  (P .  „  •  P .  .  p .  ) 


n 


j+n  "j  “A  '“j+1  j+2  ***  “j+m 

It  was  still  necessary  to  determine  if  P 


3+i 


for  i  =  1,2, • • • ,m 


was  a  P^P^,  or  matrix.  This  was  determined  as  follows: 

In  the  first  of  the  n  cycles  A  fired  his  (j+l)st,  (j+2)nd, 
L*  th 

.  ,.,(j  +  fT — )  shots  while  B  fired  his  1st,  2nd,..., 

.aB 

L*  th 

(:pr — )  shots  and  the  time  that  each  of  these  shots  occurred 

lxa 

was  determined  from  Eqs.  II-l  and  II-2. 

After  ordering  the  matrices  in  the  cycle  let  P*  =  product 

of  these  m  matrices (m-1  if  there  were  simultaneous  firings  in 

the  sequence),  i.e.,  P*  =  (P  j+l<1?  j+2*  *  *  *  *p  j+m^  *  Also  let 

A*  =  P  •  P,3.  But  since 
o  A 


Vfvt  0 


then  P.^ 
A 


1 

0 

0 


SO  that  A*  -  P0.Paj  =  (qAj  ,  PA.KA 


i-1 


,  0,  0)  —  (ai,a2,a3,a4) 

j+n,lf  pj+n,2'  pj+n,3'  pj+n,4^ * 

For- the  solution  to  the  duel,  the  individual  terms  of  Pj+n 


then  P . ,  =  A*P*n  =  (P4 

D+n 


were  de  termineu . 


The  matrix  P*  was  of  the  form 
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then 


It  was  apparent  that  the  2nd,  3rd  and  4th  terms  of  the  first 
row  of  P*n  were  geometric  summations  and  were  expressed  as: 
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Hence  Pj+n  was  expressed  as  the  following: 


Pj+n  =  A*P*n 


-  (qAj  *  oy  Pln  p2i^]  P3[^]  P4[^n-] 


n  i  1-pin  i  3  i-1  1-P 4 

(P1  qA  '  P2[l-,p1  ]  qA  +  PAi^qA  '  P3[I^[~]  qA  ' 

l-p_n 

p4 tl-pi  ^  qA3  *  * 

Since  the  objective  of  this  analysis  was  the  determination 
of  the  systems'  survival  probabilities  as  the  number  of 
transition  cycles,  n,  got  large,  the  results  were  presented 
in  the  form  of  Table  X.  As  a  quick  check,  it  was  easily 
seen  that  that  the  sum  of  the  state  probabilities  in  the 
limit  as  n-*-°°  was: 

Ai=l  A  A  l-Pi 

_  -  r 1~qA^ 1  A  „  j ,1-P1X 

PA  T-q  *  qA  1-p 


=  1  -  qAD  +  qA3  =  1,  as  it  should.  Using  the  sim- 

plified  notation  where  A*  =  (o^,  a 2,  a3r  a4)  then  Table  II 
was  derived  from  Table  I. 
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TABLE  I 


DUEL  STATE  PROBABILITIES 


Pr  (A  alive; 
B  alive) 


After  n  Transition  Cycles  As  n  +  “ 


qA3  ?in 


l  _  i-lV  P2  *A3PinP2 


Pr  (A  alive;  p  J  q  *  PA  7  < 

B  dead)  Ai=l  A  1  P1  1  P1  Ai*l 


Pr  (A  ds-ad; 
B  alive) 


1-Pl  1'Pl 


Pr  (A  dead; 
B  dead) 


TABLE  II 

DUEL  STATE  PROBABILITIES 
After  n  Transition  Cycles  As  n  +  » 


Pr  (A  alive; 
B  alive) 


Pr  (A  alive; 
B  dead) 


Pr  (A  dead; 

*1*3 

alPl  P3 

B  alive 

If 

1*pi 

Pr  (A  dead ; 

°lP4 

“lPl11  *4 

B  dead) 

i-p-T 

1-Pl 

al  p3 


B.  EQUIVALENCE  TO  PREVIOUSLY  DEVELOPED  ANALYSIS 

The  MARKOV  analysis  of  the  Stochastic  Duel  involving 
fixed  firing  rates  was  shown  to  hold  for  the  model  developed 
by  Schoderbeck  and  outlined  in  Section  I.  .In  this  model,  the 
times  between  firings,  a  and  b,  were  assumed  to  be  equal  and 
the  time  advantage  of  combatant  A  was  assumed  to  be  less  than 
a. 

Given:  a  =  b,  T<a  PQ  =  (1,  0,  0,  0) 

Then,  XA  -  XB  and  there  was  no  possibility  of  a  simultaneous 
firing  at  any  time  as  the  duel  proceeded. 
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Since  T<a,  cycle  was  of  form  DA,  so  P*  was  just  the  product 


of  the  two  matrices,  PA*PB. 


/ 


%  0  PB  0  \  /«A  Pj 


P*  = 


0  10 
0  0  1 

\  0  0  0 


0 
0 

1  / 


fop  B  pAqE  PB  °\ 


0 

0 


1 

0 

0 


0  0 


1 

0 


and 


y  n  1_(qAqB)n 

/(qAqB}  PAqB[l~qAqB  ]  PBlI^q^  ]  0 


P*n  = 


0 

0 


\  0 


1 

0 

0 


0 

1 

0 


n 


—  —  n  n  ^  ^qAqB^ 

Pl+n  ~  A*P*  =  (qAiqAqB)  '  PA  +  PAqAqB  *1-0  « 


'qAqB 


n 


lr(qAqB} 


■qAqB 


\ 


0 

1 


The  results  were  best  summarized  in  Table  III,  Checking  the 
result  obtained  by  this  analysis  and  that  developed  iri  Sec¬ 
tion  II,  it  was  seen  that  the  two  methods  were  equivalent. 

But  this  was  not  unexpected,  as  the  MARKOV  chain  analysis 
was  basically  the  conditional  probability  analysis  in  matrix 
notation. 

Analogously  the  MARKOV  analysis  was  shown  to  yield  equiv¬ 
alent  results  of  the  ANCKRR-WIL1IAMS  model  also  summarized 
in  Section  X,  as  was  expected.  The  clearest  way  to  illus¬ 
trate  this  equivalence  was  by  example,  since  the  firing  rates 
were  not  specified  in  the  ANCKER-WILLIAMS  model,  and  examples 
are  presented  in  the  next  Section. 
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TABLE  III 


DUEL  STATE  PROBABILITIES 


After  n  Transition  Cycles  as  n  °° 


Pr  (A  alive,  B  alive) 

Pr  (A  alive,  B  dead) 

Pr  [A  dead,  B  alive) 
Pr  (A  dead,  B  dead) 


n 


pa 


n 


/-'Vb’  , 

1aPb  {I^ - > 


■qAqB 

0 


p 


A 


1“qAqB 

qA*B 
1_  qAqB 


C.  EXAMPLES 

In  this  section  three  examples  are  presented  illustrating 
the  use  of  the  MARKOV  analysis.  The  first  two  examples  em¬ 
ployed  no  time  advantage  for  combatant  A  to  show  equivalence 
of  this  analysis  to  that  developed  by  Ancker  and  Williams. 

The  other  example  illustrated  the  type  of  duel  that  was  able 
to  be  analyzed  by  the  MARKOV  method  but  not  by  either  of  the 
methods  developed  by  Schoderbeck  or  Ancker  and  Williams. 

EXAMPLE  1.  PA  =  .6,  PB  =  .5,  XA  =  4  shots  per  minute 
XD  =  2  shots  per  minute,  T  =  0. 

D 


A 4  -6  0  °\ 

I  0  1 0  0  \ 

PA  =I  |] 

**  I  A  A  1  A  I 


0  0  10 
0  0  0  0 


\ 

r 5  0  .5  0  ^ 
'01  00 

r 

loo  10 

\  0  0  g  i/ 

PAB 

\ 
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(1/  0,  0,  0) 


A 


£ 

1 


>  2  W1 

A_  =  s-  =  r— 

3  1  b2 


L  =  L.C.M.  (a2,  b2)  =  L.C.M.  (1,  1)  =  1 


LX.  =  4,  LX„  =  2 

A  '  3 


L*  =  L.C.M.  (LX  ,  LXb)  =  L.C.M.  (4,  2)  =  4 


„  60  L*  _  (60) (4?  _  =  0 

"  TDT4)' (2)  30  jec‘ 


In  30  seconds  A  fired  T(g^)  =  1—L11L  =  2  shots 


and  B  fired  r  (g-j-)  =  ---3-qq--— -  =  1  shot 


i  -  t|+  n  -  tn-+  n  -  o 


"-1  -  *%rr  +  sfi  -1  -  30(15-+  §o>  -i  -  2 


L*  4 

The  number  of  P^  matrices  =  - 1)  =  —  1  =  1 


L*  4 

and  the  number  of  Pn  matrices  =  - 1)  =  —  1  =  ( 


lLX, 


and  there  was  one  P^B  matrix  in  the  cycle. 


Then , 


-  =  P.-P,3  =  P0‘PA°  =  (1,  0,  0,  0)  =  (a,,  a,,  a,  aA) 


o  A 


’I' 


3, 


P*  *  VPAB  "  -6  0  0 


0  10  0 
0  0  10 


1.2 

0 

0 


\  0  0  0  1/  \  0 


3 

1 

0 

0 


2 

0  0 
1  0 


.  3 \  /. 0  8  .72  .08  . 


\ 


0  1/  \0 


0 

1 

0 


,12\ 
0 
0 


1/ 


P .  +  n  =  P„  =■-  F  «P* 

3  no 


n 


So  as  n  +  using  results  found  in  Table  II/ 
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5  shots ,  and 


'  aA  5 

In  60  seconds  A  fired  x(g-Q-)  =  60  (g^-)  = 


-B  fired  T(g|y-)  =  60  (g~)  =  2  shots 

j  «  [1  +  i]  =  o 

>  cl 

N_1  =  t^6CT  +  6(f^  1  60  ^5o”  +  6(f*  _1  ~  6 


T  ^ 

The  number  of  P.  matrices  in  a  firing  cycle  =  (=~ - 1) 

A  LAfi 

10  ,  A  A 

=  rrmr  - 1  =  4<  ana 

t  ie  IQ 

the  number  of  Pg  matrices  =  ^ - 1)  =  (fxyjjTJ —  D'  =  ^ 


A 


and  there  was  one  Prn  matrix. 

AW 


Then, 


A*  =  P0-PAj  =  P0‘PA°  =  (1/  0,  0,  0)  =  (a^.,  a 2,  o3,  a. 


P  =  pppppp 
*  a'Wa  A  AB 


and  after  doing  this  multiplication  of  matrices, 

/ ,013  .410  .564  .013  \ 

0  10  0 

0  0  10 

\  0  0  0  1 


P*  ■ 


.and  as  n  +  »,  using  results  from  Table  II , 

a.p 


lr2 


Pr  (A  alive;  B  dead)  =  a2  +  3"n_~ 


„  (1)  (.410) 

°  +  1LV6I3 - 


=  .415 


Ancker  and  Williams  method  yielded: 


Pr  (A  alive;  B  dead)  =  .415 
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EXAMPLE  3:  E*a  =  * 8 '  PB  =  * 5 '  XA  ~  4  shots  per  minute, 


=  3  vhots  per  minute,  T  =  15  seconds. 


A 


\ 


0 

0 


8  0  0  1 

/.5  0  .5  0  \ 

10  0 

p  = 

!  o  i 

0  0 

0  10 

i 

B 

O 

o 

1  °l 

0  0  1/ 

V  0  0 

0  1/ 

.1  .4  .1  .4 


AB 


Po  =  (1,  0,  0,  0) 


XA 

4  al 

=  T  = 

XB  =  r 

^1 

°2 

L 

=  L.C.M.  (1, 

1)  =  1, 

lxa  =  4' 

lxb  = 

3 

L* 

=  L.C.M.  (4, 

3)  =  12 

T 

60  L* 

lVb 

(60) (12) 
(1)  (4) (3) 

=  60  seconds 

60 

XA 

seconds  A  fired  t  (gg-) 

-  60  <6tf> 

=  4 

shots ,  and 

xB 

B  fired  t  (g^-) 

"  60  (j|) 

=  3 

shots 

j  =  [|+1!  =  =  1 


N-1  =  +  5tr>  * 1  -  4  +  3  - 1  ■  6 


L* 


The  number  of  P*  matrices  in  a  firing  cycle  =  (j^ - 1) 


B 


12 


(1)  (3) 


-1  =  3 


and  the  number  of  Pfi  matrices  = 


,L*  i  \  _  12 

(IXA  1}  "  TTTW 


and  there  was  one  matrix.  Then, 


1 


2 


This  matrix  multiplication  yielded: 


0002 

.4888 

.1110 

.4000  \ 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

X  / 

and  as  n  -5-  «>,  again  results  were  taken  from  Table  II. 


Pr  (A  alive,  B  dead)  = 


a2  +  X- 


alE>2 


Pi 


0  L  (.2) (.4888) 
0  1-.0002 


.8980' 


D.  CHARACTERISTICS  OF  THIS  MODEL 

The  MARKOV  chain  analysis  of  the  Stochastic  Duel  was  only 
another  form  of  the  conditional  probability  models  presented 
in  Section  I.  However,  the  Markov  chain  model  enabled  the 
analyst  to  consider  duels  where  one  side  had  a  positive 
time  advantage  with  no  restrictions  on  the  length  of  the 
time  advantage.  The  model  developed  here  allowed  the  anal¬ 
ysis  of  duels  involving  different  but  fixed  rates  of  fire  by 
the  two  combatants. 

In  contrast,  Schoderbeck ' s  model  of  the  fixed  firing  rate 
duel  was  restricted  to  a  time  advantage  less  than  the  time 
period  between  two  rounds  from  one  combatant  and  further 
restricted  to  the  case  where  each  combatant  had  the  same 
firing  rate.  Ancker  and  Williams'  development  of  the  fixed 
firing  rate  duel  was  restricted  to  the  case  where  neither 
combatant  had  a  time  advantage  and  both  started  firing  simul¬ 
taneously.  Also  their  model  only  furnished  results  for  the 
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outcome  probabilities  in  the  limiting  case  after  an  infinite 
number  of  exchanges.  The  Markov  model  yielded  results  for 
the  limiting  case  and  the  case  after  a  finite  number  of 
transition  cycles,  hence  a  finite  number  of  rounds. 

Another  feature  of  this  model  was  the  ease  in  which  it 
could  be  computer  programmed.  Analysts  could  be  interested 
in  the  duel  outcome  probabilities  if  either  or  both  combat¬ 
ants  had  limited  ammunition.  For  given  firing  rates  the 
number  of  transition  cycles  before  one  combatant's  ammuni¬ 
tion  was  exhausted  could  be  computed.  Then  Pr  (A  alive? 

B  dead)  could  be  determined  at  that  point. 

Possible  extensions  of  this  analysis  along  with  uses  of 
the  fixed  firing  rate  stochastic  duel  are  presented  in  Section 


IV. 


111  *  TWO  VERSUS  ONE  DUEL -FIXED  RATES  OF  FIRE 

A.  DEVELOPMENT  OF  THE  MODEL 

In  the  model  combatant  A  dueled  with  two  other  combat¬ 
ants/  B  and  C.  The  same  assumptions  as  in  the  one-on-one 
duel  held/  but  here  A  was  assumed  to  have  two  weapon  systems 
and  fired  at  B  with  one  and  at  C  with  the  other. 

Single  shot  kill  probabilities,,  then  we  re: 

PAB  =  Pr  (A  killed  B)  (on  one  shot  at  B) 

PAC  =  Pr  (A  killed  C) 
pfiA  =  Pr  (B  killed  A) 

PCA  =  Pr  (C  killed  A) 

PBC  =  Pr  (B  killed  C)  =  0 
PCB  =  Pr  (C  killed  B)  =  0 
Firing  rates  were: 


rate 

of 

A'  S 

fire 

at 

B 

(shots 

per 

minute) 

XAC 

rate 

of 

A' s 

fire 

at 

C 

(  " 

ti 

"  ) 

XBA 

rate 

of 

B's 

fire 

at 

A 

(  " 

m 

"  ) 

XCA  = 

rate 

of 

C's 

fire 

at 

A 

(  M 

•i 

"  ) 

At  any  time  in  the  process,  then,  the  duel  was  in  one  of 

eight  possible  states: 

STATE  1:  A  alive;  B  alive;  C  alive 

STATE  2:  A  alive;  B  alive;  C  dead 

STATE  3:  A  alive;  B  dead;  C  alive 

STATE  4:  A  alive;  B  dead;  C  dead 

STATE  5:  A  dead;  B  alive;  C  alive 

STATE  6:  A  dead;  B  alive;  C  dead 

STATE  7:  A  dead;  B  dead;  C  alive 

STATE  8:  A  dead;  B  dead;  C  dead 
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Transition  matrices  associated  with  the  fifteen  possible 
occurrences  that  could  cause  a  transition  from  one  state  to 
another  are  presented  below: 


CASE  1:  A  fired  at  B 


Kb 

I: 

o 


0  Pab  0 

qAB  0  PAB 

0  10 

0  0  1 

0  0  0 

0  0  0 

0  0  0 

0  0  0 


CASE  2:  A  fired  at  C 


0 

0 

0 

0 

■1 

0 

0 

0 


0 

0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0= 


0 

0 

0 

0 

0 

0 

1 

0 
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•  CASE  3:  B  fired  at  A 


CASE  4:  C  fired  at  A 


CASE  5:  A  fired  at  B;  B  fired  at  A 


CASS  6:  A  fired  at  B;  C  fired  at  A 

^[qAB  qCA^  °  [pAB  qCi\|  °  j^AB  PCa1  0  [PAB  PCa]  0 


'AB.CA 


0 

0 

0 

0 

0 

0 

0 


dAB 
0  q 

0 

0 

0 

0 

0 


P 


CA 

0 

0 

0 

0 

0 


AB 

0 

1 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

1 

0 

0 


CA 

0 

0 

0 

1 

0 


CASE  7:  A  fired  at  C;  B  fired  at  A 


dAC  ^BA. 


-Ke 


PAC  qB^  0 


•AC- BA 


dBA 

0 

0 

0 

0 

0 

0 


dAC  ^AC 


0 

0 

0 

0 

0 


1 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 

0 


PBA 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

1 

0 


CASE  8:  A  fired  at  C;  C  fired  at  A 


f*ACqci^[PACqCa]  0  0  [qACPC^[PACPCAj 


■AC.-CA 


0 

0  10  0  0 


0 

0 

0 


0  ^ACqC^f)ACqCAl  0 


0 

0 


0 

0 


1 

0 


0 

1 


0 

0 

0 

0 


0 

0 


0 

0 

0 

0 

0 

0 


0  [qAC  P]4  [pAC  PBaJ  0  01 


0 

0 

0 

0 

0 

0 

1 1 

0 

0 


[qACPCAXE 


0 

0 


PACP| 

0 

0 


\  0  0  0 

l  0  0  0 

\  0  0  0 


0  0 

0  0 

0  o 


10  0 
0  10 
0  0  1 
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0  0.  0  0 


CASE  9:  A  fired  at  B;  A  fired  at  C 

(qAB  qAc]  [qAB  PACJ  [PAB  qAc]  [PAB  pAcj 


AB‘AC  = 


\ 


0 

0 

0 

0 

0 

0 

0 


dAB 


■AB 


0 

0 

0 

0 

0 

0 


q 


AC 


AC 


0 

0 

0 

0 

0 


1 

0 

0 

0 

0 


0  0 
0  0 


1 

0 

0 

0 


0  0  0 
0  0 
0  0 
0  0 
0  0 
1  0 
0  1 


CASE  10:  B  fired  at  A;  C  fired  at  A 

(qBA  °ooo  £p1)AgCA+qBAPCA+PBAPCAi]  o  ; 


BA  •  C  A 


0 

0 

0 

0 

0 

0 

0 


qBA  0  °o 


BA 


0 

0 

0 

.0 

0 

0 


qCA  0  0 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0 

0 

0 

1 

0 

0 


PCA  ( 


o  ( 
0  ( 
0  ( 
1  ( 
o  : 


CASE  14:  A  fired  at  C?  B  fired  at  A 


In  this  model  the  initial  state  vector  was: 

PQ  =  (.1,  0,  0,  0,  0,  0,  0,  0) 

It  was  further  assumed  here  that  A  had  time  advantages 
over  both  combatants  B  and  C  and  that  these  time  advantages 
were  equal.  In  order  to  proceed  with  the  analysis,  it  was 
then  necessary  to  compute  the  number  of  times  A  fired  at  each 
of  the  other  combatants. 

Let  T  =  time  advantage  A  had  over  B  and  C  (in  seconds) . 
Then,  the  number  of  times  A  fired  at  B  in  the  interval  -[0,  T] 
was  determined  as  follows: 
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0 


J  0 


0  0  0  0 


/ 


dAB 


l-l 


0 

1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


0  0  0  0 

0  0  0  0 
0  0  0  0 
10  0  0 
0  10  0 
0  0  10 
0  0  0  1 


Likewise  the  form  of 


was: 


0 

0 


0  0  0  0 
0  0  0  0 


0 

0 

0 

0 

0 

n 


0 

0 

0 

0 

0 


[prC 0  0  0  0 
1  0  0  0  0 

o  10  0  0 

o  0  10  0 

o  0  0  10 

o  0  0  0  1 


Ji 
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and,  then 


A*  = 


i~l  jC  32 

AB  mAC  '  mAB  ^AC.£,qAC  '  PAB  a-AQ  .  iqAB  ' 

i  x  l—l 


,  jB  jC  jB 

(q q*r  /  q*n  p 


3B  i-l  3C  i-1 

PAB  PAC.LqAB  l  qAC  '  0/  °'  0/ 

1  1  i=l 


=  (a1/a2,a3,a4,  0,  0,  0 ,  0) . 


It  followed  from  the  derivation  in  Section  II  that  the 
firing  sequence  in  the  two-versus-one  duel  was  also  periodic. 
The  combination  of  the  fifteen  possible  transition  matrices 
depended  on  the  firing  sequence,  which  was  again  dependent 
on  A*-,  AD,.,  A-,.,  and  T.  In  any  case  the  form  of  P*, 

Ad  AO  »A  OA 

the  matrix  representing  one  firing  cycle  was  of  the  form: 


p*  = 

Vi‘ 

Pj+2* 

.  .  .  *P 

j+M 

where 

j  = 

jB  + 

jC 

!  Pu 

P12 

P13 

P14 

p15 

P16 

p17 

p18  \ 

0 

P22 

0 

P24 

0 

P26 

0 

P28  j 

p*  = 

0 

0 

P33 

P34 

0 

0 

P37 

P38 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 0 

0 

0 

0 

0 

d 

1 

0  / 

V  0 

0 

0 

0 

0 

0 

0 

After  n  of  these  transition  cycles,- 

then ; 
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P16  j/ll 

n-2  j 

R.l  i 


i-1 


1 


P12P26^0  ii0P13lp223 


n 


P26  J,p22 


i-1 


i=l 


r 


+ 


0 


0 

0 

1 


0 


0 

Col.  6 


n 


Pl8i=lPl1 


i-1 


+  P12P28^0  P22 


n-2 


f  \  r  i 

j=0  i=0 


+  P13P38.An  ^AnPU  P33 


[P28i^lP221  1j 

CP38iI1P33i"'1]- 


0 

0 

0 

0 


Col.  8 
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The  probability  that  the  duel  war?  in  any  one  of  the  eight 
possible  states  after  n  transition  cycles  was  taken  from: 

Wi+n  "  K*  P*n 

=  (Pr  (STATE  1),  Pr  (STATE  2),...,  Pr  (STATE  S)) 
and  the  results  of  this  multiplication  is  presented  in  Table 

IV. 

Then,  utilizing  the  notation  for  A*,  i.e. 

A*  =  (a^,  a 2/  a^,  0,  0,  0,  0) 

the  state  probabilities  of  the  duel  were  computed  as  the 
number  of  transition  cycles  got  very  large  (n-> »  )  .  Table  V 
was  then  derived  from  Table  IV. 


TABLE  IV 


State 

Probabilities 


Pr  (STATE  1) 


Pr  (STATE  2) 


Pr  (STATE  3) 


Pr  (STATE  4) 


DUEL  STATE  PROBABILITIES 


After  n  Transition  Cycles 


jB  jC  n 
qAB  qAC  P1L 


jB  jC  9  i-1  n-i 

qAB  qAC  P12i|]pll  P22 


,  jC  i-1 

+  PAC  qAB  P22  i|1qAC 


jB~  jc  p„  I  r  i_1-  n“i 

i=l 


qAB‘'-  W  pisj  Pll"  ~P33 


3B  iej. 


+  PAB  qAC^  *  P13  ^J^AB 


n 


to3B  «*CjCtPl4iIxPu1'1 


n-2  j 


l  3-1 


P12  P24.X  ,LP11  P22 


‘3=0  i=0 
n-1  3 


+  P13  p34.^  .  ^  P11  p33^  ^ 

j=0  i=0 


jB  ^9  i-1  9  i-1 

+  pACqAB  P24i|1qAC  ^^22 


PABqAC^  P34.|nqAB  ±J1P33 


j? 


.  ,  n  •  «| 

1"1  r*  _  1 


,  ■  i-l^r  i-1 

T  PABPACi|1qAB 


46 


TABLE  IV — Continued 


TABLE  V 
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EXAMPLE: 

Let  PAB  =  *8'  PAC  =  *S'  PBA  =  *6'  PCA  =  ,5 

}.AS  =  4  shots  per  minute,  XjvC  -  2  shots  per  minute 

Ab^  =  3  shots  per  minute,  Ac^  -  2  shots  per  minute 

and  T  =  15  sec.  The  firing  sequence  is  illustrated  in  Fig. 


15  15  15  15  15 

A  at  B  —  X - X - X - -X - X _ X- 


A  at  C  — -X- 


B  at  A 


C  at  A 


20  20 
- X - X- 


•<-  T  ->■ 


Figure  3.  Firing  Sequence 


Then  P  ==  (1,  0,  0,  0,  0,  0,  0,  0) 


time 


f  j.  -  r  (4)  CL5)  .  .  , 

l^0  +  1J  -  t — 6^5 - +  1]  =  1 


A  'r 


+  1]  =  [12H15).  +  l}  m  , 


—  p  p  3^  p  jC  _  p"  p  p 

O  *AB  *AC  "  O  AB  AC 
=  (.04,  .16,  .16,  .64,  0,  0,  0,  0) 


P  •  P  •  P 
2+12+2  ’  *  '*2+6 


PAB*BA*CA  PAB*AC  PBA  ?AB*CA  PBA  PAB«AC 


Multiplying  these  matrices  yielded: 
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000001  .000045  .003896  .171393  .161062  .004147  .659456  0 


.000102 

0 

0 

0 

0 

0 

0 


388861 

.44 

1 


.131037 

0 

0 

0 

1 

0 

0 


0  .48 
.55  0 

0  0 
0  0 
0  0 
1  0 


Thus,  at  end  of  duel  (as  n->-w) 

Pr  (A  alive;  B  dead;  C  dead)  =  Pr  (STATE  4) 

alP14  ,  alp12p24p12*  ,  alP13P34P13*  ,  a2p24  ,  a3P34  , 
1"P11  (l-p12*)2  (l-p13*)2  1-p22  1~P33  4 

where  P12*  ~  raax  (Pj^/  P22)  =  *000102  and 


* 

12 

=  max  (P1]L/ 

P22 

* 

13 

=  max  (PX1/ 

P33 

Pr  (STATE  4) 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  development  in  this  thesis  of  the  fixed  firing  rate 
stochastic  duel  was  presented  to  serve  as  an  analytical 
tool  in  evaluating  present  or  proposed  weapon  systems.  The 
use  of  Markov  chains  in  the  model  enabled  the  analysis  of 
more  complex  but  also  more  realistic  weapon  system  engage¬ 
ments  than  other  models  previously  developed.  The  state 
probabilities  of  the  model  in  the  limit  as  the  number  of 
transition  cycles  increased  were  functions  only  of  single 
shot  kill  probabilities  and  hence  easily  computable  and  com¬ 
puter  programmable.  By.  computer  programming  this  model, 
parametric  studies  of  firing  rates  ,  single  shot  kill  prob¬ 
abilities  and  time  advantage  could  be  performed.  Also,  if 
the  number  of  transition  cycles  was  fixed  as  a  function  of 
one  combatant's  limited  ammunition  supply,  then  state  prob¬ 
abilities  can  be  determined  for  the  duel  where  ammunition  is 
limited. 

The  two-versus-one  duel  evolved  utilizing  the  same  tech¬ 
nique  as  the  fundamental  duel  analysis.  Results  were  more 


complex,  however  state  probabilities  still  depended  solely 


on  kill  probabilities  and  time  advantage.  It  was  thought 


that  the  two-versus-one  duel  would  be  extremely  difficult  to 


develop  as  a  conditional  probability  model,  even  without  con¬ 
sidering  a  time  advantage. 
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Stochastic  duels  involving’ combatants  with  fixed  rates’ 
of  fire  can  be  of  considerable  importance  in  the  evaluation" 
of  weapon  systems.  Using  models  developed  by  C.  J.  Ancker 
and  others,  the  fixed  rate  of  fire  duel  can  be  compared  to 
the  duel  where  time  between  firing  is  a  random  variable. 
Analysis  oC  the  random  firing  rate  and  fixed  firing  rate 
models  could  yield  the  optimal  firing  doctrine  for  given 
parameters.  To  do  this,  the  mean  values  of  the  random  fir¬ 
ing  rates  should  be  set  equal  to  the  fixed  firing  rates  in 
the  Markov  chain  model. 

B.  RECOMMENDATIONS 

Further  analysis  into  the  fixed  rate  of  fire  stochastic 
duel  can  consider  the  following  areas: 

1.  Analysis  of  the  duel  considering  properties  inherent 
to  discrete  parameter  Markov  chains  found  in  Parzen  [Ref.  4] 
and  specifically  investigating  such  properties  as: 

a.  First  passage  probabilities  and  first  passage 

times . 

b.  Absorption  probabilities  and  mean  absorption 

times . 

c.  Stationary  distributions. 

d.  Limiting  occupation  times. 

2.  Analysis  of  the  two-versus-one  duel  where  the  time 
advantage  combatant  A  has  over  B  and  that  which  A  has  over 
C  are  unequal. 

3.  The  distribution  of  rounds  fired  in  a  stochastic  duel 
utilizing  techniques  of  C.  J.  Ancker,  Jr.,  and  A.  V.  Gafarian 
[Ref.  3] . 
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APPENDIX  A 


S  QME  CON  VE  RGLXCr.  P POOPS 


A-I.  lira  ?  p  i-1  _  n-i  , 

n*'“  i=l  11  '  ?kk  =  °?  k  =  2/  3? 


Pll'  Pkk 


PROOF:  Let  pu  >  and  substitute  pn  for 
the  expression  above.  This  yields: 


lim  Y  „  n-1 

n->™  >  Pn 


Applying  L 1  Hospital's  rule: 


p-yuo  n  Pjj 


A-II , 


lira  ,  n 


f _ 11  x 

n->#0  ~-n+F  = 


-  p  -  n+1 
'1 1  1»  p 


-]  =  0 


lim  nY2  ?  „  i  j-i  Plk' 

°  A*  A  P11  P\rA  =  — — 


n*^00  .  i*  Pii 

1=0  1=0  x  KK 


(i-p,,.*) 


k  =  2'  3;  P1I'  pkk  <  1 
where  plfc*  =  max  (pu,  Pfck, 

PROOF:  Substitute  Plk*'for  pu  and 
expression.  This  yields: 


pkk  xn  akove 


lim  "f2  }  p  -  lira  nF2  •  j 

n  j=0  i=0  lk  '  n«  .1  3  Pjk*3 


n-2 


.Lim 

n+®  ^lk  dpjk*  '^lk 


P,v*  l 


T  (Pn.*3) 


,  n"2  4 

d  r  -  *D 


lim  p  *  ...JL-..  y  p 

n->®  pik  dp^k*  1  lk 


i-Piv* 


n-1 


li»  p  *  f Hite 

n*>»  *■  ik  dpkk*  1  l-p^f 


lim  *  (1-Plk*)  t  -  (n-l)p  k*n"2]-(l-p  *n-1)  (-1) 

n-yco  Pik  [ - ±iS - ±h - - - iiS - ] 


(1  “  Plk*} 


2 


lim  r  ^lk 
n+®  ,,  ..  ^2 


Pik* 


a-Plk*)‘  u-plk*> 
p  * 


2  Un^lp^*"-1) 


lk  ?  Kn-DPlk*n-2]) 


(i-Pj.,*) 


(A-II-1) 


but  the  second  and  third  terms  of  Eq.  A-II-1  went  to  zero  in 
the  limit  as  n->®  by  the  proof  of  A-I  above. 
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